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We analyze the coupling between the internal degrees of freedom of neutron stars in a close binary, 
and the stars’ orbital motion. Our analysis is based on the method of matched asymptotic expan¬ 
sions and is valid to all orders in the strength of internal gravity in each star, but is perturbative 
in the “tidal expansion parameter” (stellar radius)/(orbital separation). At first order in the tidal 
expansion parameter, we show that the internal structure of ea ch star is unaf fected by its compan¬ 
ion, in agreement with post-l-Newtonian results of Wiseman (5r-qc/970401E). We also show that 


relativistic interactions that scale as higher powers of the tidal expansion parameter produce quali¬ 
tatively similar effects to their Newtonian counterparts: there are corrections to the Newtonian tidal 
distortion of each star, both of which occur at third order in the tidal expansion parameter, and 
there are corrections to the Newtonian decrease in central density of each star (Newtonian “tidal 
stabilization”), both of which are sixth order in the tidal expansion parameter. There are additional 
interactions with no Newtonian analogs, but these do not change the central density of each star up 
to sixth order in the tidal expansion parameter. These results, in combination with previous analy¬ 
ses of Newtonian tidal interactions, indicate that (i) there are no large general-relativistic crushing 
forces that could cause the stars to collapse to black holes prior to the dynamical orbital instabil¬ 
ity, and (ii) the conventional wisdom with respect to coalescing binary neutron stars as sources of 
gravitational-wave bursts is correct: namely, the finite-stellar-size corrections to the gravitational 
waveform will be unimportant for the purpose of detecting the coalescences. 


04.25.-g, 04.40.Dg, 97.80.-d, 97.60.J 


I. INTRODUCTION AND SUMMARY 

Recent numerical simulations by Wilson, Mathews and 
Marronetti of the late stages of inspiral of neutron star 
binaries have predicted the following surprising result: 
The individual neutron stars are apparently subject to 
a crushing force of general-relativistic origin which can 
cause the stars to collapse and form black holes, before 
they reach the dynamical orbital instability that marks 
the end of the inspiral These numerical simula¬ 

tions were fully relativistic, but assumed a conformally 
flat spatial metric, and also employed an approximation 
scheme in which the gravitational field was constrained 
to be time-symmetric at each time-step in the compu¬ 
tation. These approximations and assumptions give cor¬ 
rect results for spherically symmetric systems and also to 
the first post-Newtonian approximation |^; beyond this, 
however, their domain of validity is not well understood. 

The Wilson-Mathews-Marronetti prediction is in dis¬ 
agreement with other, independent, fully relativistic, nu¬ 
merical simulations which employ similar approximations 
iQ,^, with post-l-Newtonian numerical simulations [^, 
and with Post-Newtonian 0-i and perturbation 0 cal¬ 
culations, which we discuss further below. Therefore it 
seems likely that star crushing does not occur in reality, 
although the issue is still somewhat controversial. 

The star-crushing scenario, if correct, would have pro¬ 
found implications for the efforts to detect gravitational 
waves produced by neutron star inspirals with gro und 
based interferometers such as LIGO and VIRGO [|ll| . A 


crushing force which is strong enough to cause an insta¬ 
bility to radial collapse of a neutron star would consti¬ 
tute a strong coupling between the orbital motion and 
the internal modes of each star ||^ , and would transfer a 
substantial amount of energy from the orbital motion to 
each star. Specifically, let ^instability be the value of the 
orbital separation L at the onset of instability, let R be 
the initial radius of the neutron star, and let AR be the 
amount by which the star’s radius is decreased before the 
instability occurs. Then the ratio of the energy AEgtar 
absorbed by the neutron star ||^ to the energy AEgw ra¬ 
diated in gravitational waves between L = Tinstabiiity/2 
and L = ^instability would be approximately 

AEgtar / nstability \ fAR\ 

AEgw ^ V R )\^J' ■ 

The Wilson-Mathews-Marronetti simulations predict 
that an instability occurs at ^instability ^ 5 R, and that 
AR/R ~ 1/20 Q. Therefore the ratio ( |l.lD is predicted 
to be of order unity, and the crushing effect gives rise to 
an order unity perturbation to the inspiral rate of the 
orbit and to the phase evolution of the emitted gravita¬ 
tional waves 1^. Since high accuracy theoretical tem¬ 
plates are required in order to extract the signal from 
detector noise, the star-crushing scenario would imply 
that currently envisaged search templates (which are 
calculated neglecting all orbital-motion—internal-mode 
couplings) would need to be completely revised. There¬ 
fore, it is important to find out whether or not the star 
crushing effect occurs. 
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In the Newtonian approximation, the coupling between 
the stars’ orbital motions and their internal motions has 
been analyzed in detail BH- The Newtonian coupling 
is very weak, too weak to affect the gravitational wave 
signal except in the last few orbits before coalescence 
|p7|-^. Moreover, the Newtonian interaction energy is 
2 — 3 or ders of magnitude smaller at L ^ 5R than the 
amount (0) which would be required to crush neutron 
stars when the necessary AR/R is of order several per¬ 
cent Although post-Newtonian couplings have not 

yet been analyzed in detail Q , the expectation has gen¬ 
erally been that post-Newtonian or relativistic couplings 
would simply be of the form (Newtonian coupling) x 
0(M/R) or 0{M/L), where M is the neutron star mass. 
That is, only small fractional corrections to existing New¬ 
tonian couplings are expected. However, the Wilson- 
Mathews-Marronetti prediction suggested instead that 
relativistic couplings could dominate over the Newtonian 
ones, and highlighted the need to understand the details 
of these couplings. The purpose of this paper is to explore 
and elucidate the relativistic, post-Newtonian couplings. 

In the simulations of Wilson, Matthews and Mar- 
ronetti Q , the central density pc of each star was seen to 
increase during the inspiral, and the instability to radial 
collapse occurred once the fractional increase in central 
density reached a few percent. In this paper, we shall 
not consider the issue of stability to radial collapse, but 
instead, following Refs. [^|Jl^], we shall focus attention 
on how the central density changes during the inspiral. 
Focusing on the central density allows us to investigate 
the existence of radial crushing forces (which presumably 
cause the instability to collapse in the simulations). 


A. Relativistic coupling between orbital motion and 
internal degrees of freedom in neutron star binaries 


For a neutron star in a binary system, let L, i? and M 
be the orbital separation, stellar radius and stellar mass, 
respectively. Then, there are two natural, independent, 
dimensionless parameters that characterize the system: 
the strength of internal gravity in each star 


M 



( 1 . 2 ) 


which is of order ^ 0.2, and the tidal expansion param¬ 
eter 


R 



(1.3) 


which gradually increases during the inspiral. In terms of 
these parameters, the post-Newtonian expansion param¬ 
eter for the orbital motion is Corbit = M/L = ea. The 
standard post-Newtonian approximation scheme consists 
of expanding in e and Corbit, treating these quantities as 
formally of the same order. In the context of a neutron 
star binary, the usual terminology for describing the size 


of gravitational effects (Newtonian, Post-l-Newtonian, 
Post-2-Newtonian etc.) is somewhat ambiguous, since 
an effect of post-n-Newtonian order could scale as e“(ae)*' 
for any a, b with a -I- 5 = n. In addition, some authors 
would classify an interaction which scales like as be¬ 
ing of effective post-Newtonian order t, irrespective of 
the value of s, since the index t controls the perturba¬ 
tion to the phase evolution of the emitted gravitational 
waves. Therefore, in this paper we will not use the post- 
Newtonian terminology. 

Let us start by reviewing the situation in Newtonian 
gravity, at order (see Table | below). As is well known, 
each star distorts the other at 0{a^) via the quadrupole 
tidal interaction, and the second order response of each 
star to the quadrupole tidal field of its companion gives 
rise to a decrease in its central density pc at 0(a®) [Q. 
Each star expands slightly and thus becomes more stable 
against radial collapse 0. 

Consider now the situation at higher orders in e. If 
there were a relativistic crushing force that changed the 
central density of each star, then this change in density 
6Pc would scale in some way with the dimensionless pa¬ 
rameters a and e: 


Spc 

Pc 


oc a* e'*, 


(1.4) 


for some integers or half-integers s,t > 0. Such an ef¬ 
fect would dominate over the Newtonian effect at large 
orbital separations if < < 6. Note that naive arguments 
(which turn out to be incorrect) do suggest scalings of the 
form (0) with low powers of a [Q. In the numerical 
simulations it was seen that 6pc/pc oc |Q, while the 
scaling with e was not clear. A post-I-Newtonian calcu¬ 
lation by Wiseman [0 has shown that there is no change 
in Pc at the order 0(a e ), but did not rule out the pos¬ 
sibilities 5Pc/Pc oc ea* with 1 < t < 6, or SpdPc oc e® a 
with s > 1. 

A different approximation scheme was used by Brady 
and Hughes to investigate relativistic interactions [ p^ . 
Suppose that the two stars are labeled A and B, and that 
their masses and radii are Ma, Mb, Ra and Rb- Brady 
and Hughes showed using perturbation theory that there 
is no change in the central density of star A to linear 
order iw Mb- However, this result would not rule out a 
change in central density scaling as, for example, 


bpc MaM"^ 

Pc RaL"^ 


(1.5) 


If we specialize to the equal mass case Ma = Mb and 
Ra = Rb, the scaling (1^) reduces to Spc/pc oc ea^. 
Thus, the Brady-Hughes analysis does not necessarily 
rule out the type of behavior seen in the numerical sim¬ 
ulations. 

In this paper we analyze the change in central density 
to all orders in e, but perturbatively in a. We show that 
in a binary system, the excitation of the internal degrees 
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of freedom of a fully relativistic spherical star is quali¬ 
tatively the same as that of a Newtonian star. At order 
0{a^), the star responds linearly to the external gravi¬ 
tational field, and is distorted; but there is no excitation 
of the stars spherically symmetric, radial modes. At or¬ 
der 0 (a®), the second-order response of the star to the 
external field generates an excitation of the stars spher¬ 
ically symmetric modes and a corresponding change in 
central density. We deduce that to all orders in e, any 
changes in central density must scale in the same way as 
in Newtonian theory: 


^Pc 


( 1 . 6 ) 


as a —> 0; see Table | below. Our result is restricted 
by the assumption that the neutron stars are not spin¬ 
ning. This restriction is unimportant since in Ref. 0 the 
crushing effect is seen when the neutron stars have very 
small net spins. Note that the result (O) is not based 
on an analysis that includes “only tidal interactions”. 
Instead, our analysis is a demonstration that there are 
no general-relativistic interactions that contribute to the 
leading-order change in central density other than the 
Newtonian-ty pe t idal interactions. 

Our result (1.6) is in disagreement with Refs, 
the regime a 1 where there should be agreement: Fig. 
1 of Ref. 1^ shows that the scaling found in the numerical 
simulations is SpdPc oc as a —> 0. This disagreement 
in scaling strongly suggests that the star crushing effect 
seen in Refs. |]^,|| is not physical. 

Our derivation of the result (a takes place in two 
stages. First, in Secs. i 0 and 0 we analyze the 
change in central density for a star moving in a fixed, 
external, vacuum gravitational field. We show that in 
this context 5pc scales as 7?.“^, where TZ is t he r adius of 
curvature of the external spacetime [cf. Eq. (4^) below]. 
In Sec. ^ we extend the analysis to two stars in a binary, 
and deduce the scaling (|l. 6 | ). 

We use units in which the speed of light c and Newton’s 
gravitational constant G are unity, and use the sign con¬ 
ventions of Ref. I^j. Indices a,b,c,... will be abstract 
spacetime indices in the sense of Wald [Q, thus equa¬ 
tions involving such indices will be valid in all coordinate 
systems. Indices a, (3, 7 ,... and i,j, k,... will be con¬ 
ventional indices, the former running over 0,1,2,3, the 
latter running over 1,2,3. 


II. INTERACTIONS OF A FREELY FALLING 
BODY WITH AN EXTERNAL GRAVITATIONAL 
FIELD 


In this section and in the following two sections we shall 
consider a neutron star moving in some arbitrary back¬ 
ground vacuum gravitational field, for example a super- 
massive black hole. The key technical tool in our analysis 
is the method of matched asymptotic expansions, as ex¬ 
plained in, for example, Refs. ]2^-|^. This method has 


been used in the past in general relativity primarily to 
derive equations of motion for bodies moving in external 
gravitational fields |^,^. However, the method also 
lends itself naturally to analyzing the effect of an exter¬ 
nal gravitational field on the internal structure of a fully 
relativistic, self-gravitating body. The key feature of the 
matched asymptotic expansion method is a separation 
of lengthscales/timescales: the radius of curvature of the 
external spacetime is assumed to be much larger than 
the lengthscales characterizing the neutron star, and the 
timescales over which the external curvature is changing 
(as perceived on the neutron star’s worldline) is much 
longer than the internal dynamical timescale of the neu¬ 
tron star. 

The methods and result which we discuss in this sec¬ 
tion and in Secs. Ill and 0 below do not apply directly 
to two stars of comparable mass in a binary, since in that 
context there is no external, fixed, background gravita¬ 
tional field. In Sec. ^ we show how to mesh the matched 
asymptotic expansion method with the post-Newtonian 
expansion method, and extend our analysis to be appli¬ 
cable to two stars in a binary. The discussion of motion 
in fixed external gravitational field of Secs. i -0 is in¬ 
cluded as background and motivation for the analysis of 
Sec. 0. 


A. Constructing the spacetime: setup 


Let the metric of the background field (for example, 
the supermassive black hole) be g^\ and suppose that 
there are no matter sources in the region of interest in 
the background spacetime, i.e., it is a solution of the 
vacuum Einstein equations. Suppose also that F is some 
geodesic in this spacetime. Then, one can pick Fermi- 
normal coordinates = {t, x^) adapted to this geodesic, 
such that F is the curve a:® = 0, and such that the metric 
is of the form 

9ab — Vab + ("^- 1 ) 

where fjab is the flat, Minkowski metric fjab = 
— (dt)a{dT)i,+5ij{dx3')a(dxdb, and where is quadratic 
in distance from the geodesic. More specifically, near the 
geodesic F we have IM 


h^apdx°‘dx^ 


-Roiomxd"^ dP - ^Roijmx’’x'^ dtdx^ 
-^Riijmxd'^ dx'^dx^ -l- 0 (|i®p), ( 2 . 2 ) 


where the various Riemann tensor components are eval¬ 
uated at {t,xd = i.e., on F. 

We next introduce some notation. Let TZc denote 
the order of magnitude of the radii of curvature of the 
background spacetime along F, so that TZ^ = typical 
value of the Riemann tensor components on the right 
hand side of Eq. (2.2). Let C and T be the lengthscale 
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and timescale, respectively, over which the curvature is 
changing, given schematically by Rajj-fijC, ~ ^iRap-yS 
and Rajs^sjT ~ '^oRap-fS- Below we will treat TZ^, ^ 
and T as formally of the same magnitude, and will de¬ 
note these lengthscales collectively by TZ. 

Consider now a completely different spacetime: a 
static, spherical, isolated neutron star, which we model 
as a perfect fluid obeying some equation of state p = p{p) 
Let its mass be M and its Schwarschild radius be R. 
In a suitable coordinate system a;“ = (t,a;*), the neutron 
star metric can be written as 


„(NS) _ „ , ;,(NS) 

iJab — 'lab + 'I'ab > 


(2.3) 


where pab = -{dt)a{dt)b + 5ij{dx'')a{dx^)b, r = y/5ijX'x^, 
and at large r 

h%^Ux°‘dx^ = - [-2M/r + 2M^/r^] dt^ 


2M 3M2 


+0 


r 2r^ 

1 


5jkdx'’dx^ 


(2.4) 


T he t ask n ow i s to construct, starting from the met¬ 
rics & and (|2.3|) , an approximate solution of Einstein’s 
equations that represents the neutron star traveling along 
the curve T in the background spacetime, in the limit 
where R TZ. In Sec. IIB we describe, without proof, 


the resulting spacetime to lowest order in the dimension¬ 
less parameters 


and 


7 i = R/TZ 


72 = M/TZ. 


(2.5) 


( 2 . 6 ) 


This leadin g ord er spacetime is well known. Then, start¬ 
ing in Sec. 11C we describe a systematic procedure for 
calculating the spacetime metric to successive orders in 
the parameters (^) and ( ^) , which can be used 

to verify the results of Sec. P^IB| )^ Finally in Sec. 0 we 
deduce the scaling of the change in central density of the 
star. 


B. The leading order spacetime 

Consider the region R r TZ in the background 
spacetime (2.2), which we will call the matching region. 
In this region, consider the metric 


(HI) is an approximate solution of the vacuum Einstein 
equation in the matching region, since and 

are both approximate solutions of the linearized vacuum 
Einstein equation. To l ead ing or der in the dimens ionl ess 
expansion parameters (2.5) and (2.6), the metric (H) is 
the correct, physical metric in the matching region. 

What of outside the matching region? The physical 
metric of the spacetime, again to leading order in the 
tidal expansion parameters, can be obtained as follows. 
In the region r ~ i?, which we will call the interior region, 
the metric is 


(INTERIOR) 

9ab 


— Pab + 


"ab + "ab ' 


( 2 . 8 ) 


Here, the tensor is the full, nonlinear metric pertur- 

XK / dn 

bation of the neutron star. The tensor is a linearized 
metric perturbation which describes the leading order ef¬ 
fect of the external tidal field. This metric perturbation, 
together with some linearized perturbation to the star’s 
fluid variables, is a solution of the coupled Einstein plus 
perfect fluid equations, linearized about the neutron star 
background, and with all time derivative terms dropped 
[^ . [The time derivative terms scale as 1/T ~ 1/7^ 
and thus e nter at higher ord er in the small parameters 
(|2.5|) and (|2.6| ); see Sec. 11C below and Ref. @].] The 
solution is chosen such that, at large r, ^ the 
quadratic expression on the right hand side of Eq. (^.2|). 
These perturbations (of the metric and of the fluid vari¬ 
ables) describe normal mode deformations of the star 
responding adiabatically to the external tidal field. The 
boundary condition at large r determines both the met- 
ric perturbation and the perturbations to the fluid 
variables. The situation is analogous to that in Newto¬ 
nian gravity as analyzed in detail in Ref. [p0| , where the 
boundary condition on the Newtonian potential at large r 
determines the solution of the stellar perturbation equa¬ 
tions. 

In a similar way, the spacetime metric in the exterior 
region r ^ TZ can be written as 


(EXTERIOR) _ , £(NS) 

9ab - dab + 


h 


(B) 

ab ■ 


(2.9) 


Here, is of order unity, and is a linearized solu¬ 
tion of the vacuum Einstein equations on the background 
9ab ^ exterior region, which matches smoothly onto 

in the matching region. To leading order, is 

just the solution of the linearized vacuum Einstein equa- 
tions on the background whose source is a delta 
function along the worldline T. 


(MATCHING) , ,(NS) 

9ab = dab + 


JB) 

^ab 


(2.7) 


obtained by identifying the coordinates used in Eqs. ( |2.2| ) 
and (2.4) and by simply adding the metric perturbations. 
The superscript (MATCHING) indicates the region of 
spacetime in which the expression is valid. The metric 


C. Constructing the spacetime: general method 

In this section we describe a general scheme for calcu¬ 
lating the spacet ime metric perturbatively in the param¬ 
eters (^) and (^), based on the treatment in Thorne 
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and Hartle and in Mino et. a l. 1^ . T his general 
scheme justifies the claimed forms (|2.8| ) and (|2.9|) above 
of the leading order metrics. It also enables us to go to 
higher order, which is necessary since (as in Newtonian 
and post-Newtonian theory) it turns out that changes in 
the central density of the star are quadratic in the leading 

^ / D\ 

order tidal field . Therefore, in order to determine 
the leading order change in central density of the star, it 
is necessary to consider higher order perturbations which 
scale as 

There are three elements to the general procedure: an 
internal scheme, an external scheme, and a matching of 
the two schemes 


1. The internal scheme 


Let denote the manifold and back¬ 

ground metric ( |2.3D of the neutron star. The physical 
metric can be written as the background metric plus a 
sequence of pe rtur bations of various orders, in a general¬ 
ization of Eq. (2.8): 


(INTERIOR) _ (NS) 

3ab 9ab 

-L£4)j(4) 

"'ab 


„£(!) 1 o2/,(2) I p3t(3) 

-he -he 


0{e% 


( 2 . 10 ) 


Here is what we called above, and below we 
;(i) 


— u. The quantity e is a formal ex¬ 


will show that = 0. 
pansion parameter which can be set to one at the end 


^*2 scales like TZ * (but 
has no definite scaling with respect to M). We shall be 


of the calculation; each term scales like TZ 


working to order 0(e'*). The expansion ( |2.10[) constitutes 
a general-relativistic generalization of the usual tidal ex¬ 
pansion of stellar interactions in Newtonian gravity. In a 
similar way we can expand the stress-energy tensor Tab 
of the neutron star as 


rp _ rpiO) . rpy-i-J , , ^^rpy-tf 

-rab — -bnh “T I -|- fc J- I -t t -L„f, -|- fc 1 ab "T . . . . 


dl) 


( 2 ) 


Jirp{T) . 4rp{4:) 


^ ab 


^ ab 


^ ab 


ab 


^ ab 


( 2 . 11 ) 


Here is the stress-tensor of the static, spherical, un¬ 
perturbed neutron star, and the higher order terms £*T^^ 
scale as TZ~^. 

The terms with t = 1 in Eq. (2.1C) and ( ^.11 ) actually 
vanish identically. This is because when one carries out 
the matching procedure described below to determine the 
solutions, one finds that there are no pieces of the inter¬ 
nal solution that scale as TZ~^ |26|. Henceforth we shall 


anticipate this result and set = 0. 

Now the perturbations to the interior metric and to 
the neutron star are driven by the external gravitational 
fields which vary over a time scale T ~ 7?. oc 1/e. This 
fact needs to be built into the approximation scheme we 
use to derive the perturbation equations of motion. We 
can illustrate the nature of the required approximation 


scheme with the following example. Consider a scalar 
field d) on flat spacetime obeying the equation 


□*=l-^ + v 


d>(x,t) = p(x,t). 


( 2 . 12 ) 


where the source /9(x, t) is of the form p(x, t) = et), 
and the function po is independent of e. If we define 
$o(x,t) = <i)(x,t/e) and r = et, Eq. ( ^.12 ) can be re¬ 
written as 




$o(x,r) = po(x,r). 


(2.13) 


If po is now specified as a power series in e, we can solve 
for <I>o order by order in e using Eq. (Hi- The equa¬ 
tion of motion (2.13) can also be obtained (dropping the 
subscri pts 0) simply by multiplying each time derivative 
in Eq. (M by e. 

In a similar way, the equations of motion of the internal 
scheme can be obtained by writing out the Einstein and 
fluid equations in terms of the contravariant components 
gab and Tab of the metric and stress tensor, substituting 
in the expansions (2.10) and (2.11), and by multiplying 
each time derivative by e. More specifically, for each e, 
there will exist an e-dependent coordinate system (i®, t) 
which varies smoothly in e and which coincides with the 
coordinate system {x^,t) of Eq. (2.4) at e = 0, such that 


,i)dx°‘dx^ = h^*g{x^, ei)dx°'dx^ (2.14) 




and 


T^^{x'‘,i)dx°‘dx^ = T^^{x\ei)dx°'dx^ (2.15) 


for t = 1, 2,3 .... Here the tensors and T//® are in¬ 
dependent of e. The Ei nstein and p er fect fluid equa - 
tions combined with Eqs. ( ^.10 ), (2.11), ( ^.14 ) and ( 2T5| ) 




now yield a system of equations for the tensors and 


analogous to Eq. ( 2.13| ). These equations will only 
be valid in the res trict ed clas s of c oordinate systems for 
which the ansatz ( ^.14| ) and ( ^.15 ) are valid. 

The approximation scheme can also be described in 
coordinate invariant terms in the following way. Let t(e) 
be a scalar field and t“(e) a vector field with 


= 1 , 


(2.16) 


whic h reduce to the t and d/dt of the coordinate system 
(^.4D at e = 0. Define the tensor 


Kie)=5^a+[]-l]t'’Vat 


(2.17) 


Let Xe '■ be the mapping which moves 

any point Ine units along integral curves of the vector 
field i.e., in suitable coordinate systems {x^,i), Xe 
maps (i®, i) to (a;®, ei). We make the following ansatz for 
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the form of the metric stress tensor 

Tab- 

^(INTERIOR) = Kmt{e)XeM (2.18) 


where and scale as T p, i.e., contain p time 

derivatives, for p = 0,1,2. We write the derivative op¬ 
erator for the metric -I- £hab, for any perturbation 


and 


where 


Tab{e) = K{e)kt{e)xeXd. 


gab = + <^Kb + s'" Kb + 

+eK2+0{e^) 


( 1 ) 


,2r(2) . ^3 l(3) 


and 


Tab = T® + eTll> + + 0(^")- 

( 2 . 21 ) 

Here the tensors Kl and are independent of e and 


s(i) 


=.(2) 


3^(3) I 4/^(4) 


(2.19) 


( 2 . 20 ) 


^ah 


^ ab 


Xe* is the pullbac k ma p. Equat ions ( ^.18 ) - ( |2.21) to - 
gether with Eqs. ( 2.10| ) - ( 2.1l| ) reduce to Eqs. ( 2.14 ) 
and ( p.l5 ) in suitable co ordin ate sy stem s. [Note that the 
transformation in Eqs. (2.18) and ( ^.19| ) leave the back¬ 
ground quantities and invariant.] Finally, the 

Einstein and perfect fluid equations 


ya ^ y(0)a -k + OK)- (2.27) 


[h,h] 


This equation defines the operators Vj^j“ and 
which are multiplicative and not differential operators 
despite the notation. Finally, we use the expansion 


y(t) a _ y(t,p) ( 

p—0 


(2.28) 


for t = 0,1,2, where contains p time derivatives 

(either acting on h or as differential operators) and thus 
scales as T~p. 


Using these notations and the expansions (2.10) and 
o, we obtain from Eq. ( |2.23D the perturbed fluid 
equations 


Y7(0,0) arpi‘2) —(1,0) o—(0) p. 

-'ob [)j(2)] ^ab — 


YjiOyO) arp{3) . ^(1,0) ar^(O) _ -r-iS) 

^ ^ab 'S' ''[^(3)] S-ab — Sb ’ 


(2.29) 


(2.30) 


Gobb, 


(INTERIOR)i 


cd 


= 0 


v“r„6 = 0, 


( 2 . 22 ) 

(2.23) 


when expanded order by order in e, yield a system of 
elliptic equations for the tensors and T^^f ■ 

Below, we shall for simplicity set e = 1 after having 
derived the equations of motion. When e = 1 we have 


and = Ky, and thus we can write the 
equations in terms of the fields and . 

In order to write out explicitly the resulting equations 
of motion, we introduce the following notations. We 

(NS) 




n(t) 


t(*) 


ah 


define the operators G^^^ [/i; and G^^^ h,h;g 

which act on metric perturbations and pairs of metric 
perturbations respectively, to be linear and quadratic 

parts of the Einstein tensor of the perturbed spacetime 
^(NS) _|_ 


Gabb^^®) +e/i] = Gab[9^^^^]+£Giy[h;g^^^^] 
+e^G^^J[h,h-,g<^^^^]+OK). 


(2.24) 


[Here and below we drop some of the tensorial indices for 
ease of notation]. We also use the expansions 


and 


y ( 0 , 0 ) arpG) 


— (1,0) a —(0) .pC4) 

■ Vj^(4)] J-ab - -^b > 


where the 4-force densities and 


are 


'r*(^) _ _^ 7 ( 14 ) <3,rT-i(0) _ ^ 7 ( 04 )^ 17 ^^^) 

^ ^ o^b ^ ab 


and 


(2.31) 


(2.32) 


^(4) _ 

■^b — 


_a^( 0 ) _ ^ 7 ( 14 ) a.Tn( 0 ) _ ■^( 0 , 1 ) arp{^) 

[^(2)] ah [^(3)] ah ^ ^ ah 




rpig) Yy(4,0) ^rp{‘2) /n 00\ 

[/t(2),b(2)]^ab '^[b(2)] ab ■ ( -OO) 


Note that the actual fluid equations of motion are ob¬ 
tained by using the formula Tab = {p + p)uaUb +pgab, by 
assuming expansions of the form ( p.ll ) for the den sity p 
an d fou r velocity u“, and by substituting i nto E qs. ( 2.2E ) 
- ( 2.31 ). However, the schematic form ( 2.29|) - ( 2.31 ) 
will be sufficient for our purposes. 

Similarly, from the Einstein equations we obtain the 
equations of motion for the metric perturbations for 


TO = 2, 3,4 : 


[h] = Giyf [h] + GiY^ [h] + GiY^ [h] (2.25) 


G' 


(1.0) 

ab 




= 8ttT 


( 2 ) 


ab ’ 


(2.34) 


and 


G^ab [h, h] = Gil’°) [h, h] + Gif) [h, h] + Gif) [h, h ], 


(2.26) 


G 


(1,0) 

ab 




= 8^T’4-Gif) 


hKg^ 


(2.35) 
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and 


G 


(1.0) 

ab 




_^(2.0) 

^ab 









solution of the Einstein plus perfect fluid equations, lin¬ 
earized about the static neutron star background, with all 
time deriva ti ve te r ms dr opped. The higher order equa¬ 
tions ( ^.2S| ), ( ^.301 ), ( ^.35| ) and ( 
structure, but contain additional source terms. We shall 
be interested in very general solutions of the equations 
which are allowed to diverge as r ^ oo; the physical so¬ 
lutions will be determined by matching to the metric of 
the external scheme. 


2.36) have the same basic 


2. The external scheme 


Let denote the manifold and metric dH) 

of the external spacetime. The physical metric can be 
written as this background metric plus a sequence of 
pert urba tions of various orders, in a generalization of 
Eq. (p[): 


(EXTERIOR) _ (B) 

9ab 3 ab 


+ eh''2+e^h''^^+0{e^). (2.37) 


Here as above e is a formal expansion parameter which 
can be set to one at the end of the calculation; in 
Eq. ( 2.37 ), each term scales like M® but has no 


definite scaling with respect to TZ. [Here we treat the 
quantities M and R as formally of the same magnitude, 
so that the expansion parameters ( |2.5D and (2.6) coin¬ 
cide 


(NS) 

ab 


^]. The perturbation was denoted above as 
The equations of motion satisfied by the metric 


perturbations and are obtained from the vacuum 
Einstein equation: 


^ab 


Gll’ [S“>;9S' 


= 0 
= 0 . 


(2.38) 

(2.39) 


The solutions to these equations will diverge as r —> 0, 
and will be determined by matching onto the internal 
solutions. 


3. The matching scheme 


In t he matching region i? <C r 77., the interior metric 
(^.10) can be written as the double expansion Q: 


(INTERIOR) , f(s.t) ( ^\ f f V 

9ab =3ab + }^}^K,> l—j . 

s=0 t=0 ^ ^ 

(2.40) 

Here the terms with t = 0 are the expansion in powers 
of M/r of the neutron star metric perturbation 
and the terms with t = 1 are the expansion in pow¬ 
ers of M/r of the metric perturbation ■ The terms 
with t = 2 give the expansion in powers of M/r of the 


metric = 


-AB) 


, which describes the leading order 


tidal gravit ation al field, and so forth. The form of the 
expansion ( ^.40 ) can be regarded as an ansatz which is 
validated by explicit matching calculations to the first 
few orders in s and t |2^; however, at higher orders 
in s and t it might be necessary to include powers of 
{M/r)log[M/r] or (r/77) log[r/77.] |Q. Moreover, the 
expansion is merely asymptoti c and is not expected to 
be convergent |^. Erom Eqs. ( p.34 ) - ( p.36 ), each term 
/rY{r/TZy in Eq. ( 2.401 ) obeys a linear elliptic 
equation with nonlinear source terms generated by the 
terms * \M/rY (r/77)‘ with s' < s and/or t' < t. 
The solutions for t > 1 diverge as r ^ oo. 

In a similar way, we can perform a double expansion 
of the external metric (2.37): 


(EXTERIOR) - , 7 

9ab =9ab + 2_^2^h 

s=0 t—0 


(s.t) 

ab 




{£)'■ 

(2.41) 


where = dijX'xT Here the terms of a given order s 
give the expansion in powers of f/77 of the metric pertur¬ 


bation h 


, while the terms with s = 0 are the expansion 

/ D\ 

of the background metric . 

The basic ide a now is t o dem and consistency between 
the expansions ( 2.40| ) and (2.41). Before this can be done, 
however, one must specify an embedding 




(2.42) 


of the neutron star spacetime into the external spacetime, 
i.e., a mapping between the asymptotically Lorentzian 
coordinates x°‘ of and the Eermi-normal coordi¬ 

nates x°‘ of M^^'>. Let us write this mapping as 

5“ =x°'{t,x') 

= 11 (D'EE "■'■’nwa.'', (2.43) 

P^O ^ / r^o I-r 

where R denotes the m ulti-i ndex (Ii,..., v) and = 
x''^ ... a;*’". Equation ( ^.43| ) is a Taylor expansion of p 
in terms both of the spatial coordinates a:* at each fixed 
t, and also in terms of the parameter M/77.. Thus, the 
worldline of the center of the neutron star gets mapped 
onto the worldline 
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X' 






The matching procedure described below can be used 
to show that the first term in Eq. (2.44) represents a 
geodesic in the background metric ', and that the sec¬ 
ond term is the first order correction to geodesic motion 
due to radiation reactio n . 

Using the expansions (|2.41|) and ( 2.43| ) we can calculate 
the pullback the external metric to the 

manifold and expand it as: 


(EXTERIOR) , 

V’*ffab =Vab + 2 ^ 2 ^ Kl 

s=0 t—0 


f) 


/r^\ t 

V7^/ 


(2.45) 


Here each term h 


/{s,t) 
ab 


can depend on the quantities 


^ 0 < s' < s and 0 < t' < t, and also on the 


embedding functions with 0 < r < t -|- 1 and 

0 < p < s. For al l choices of these embedding functions, 
the metric ( 2.45 ) will satisfy the perturbative vacuum 
Einstein equation, since it is the pullback of a metric 
which does so. 

We can now describe, schematically, the matching pro¬ 
cedure. First, solve for the general solutions for the 
metric pert urba tions u p to a given order m and I in 
both Eqs. ( 2.40| ) and ( |2.41 ). Second, specify the em¬ 
bedding ( 2.43D up to the required order in r and p, and 
calculate the expansion coefficients in Eq. ( ^.45 ). 

Third, demand that the expansions ( 2.4C| ) and ( 2.45| ) 
agree, which determines both the metric perturbation so¬ 
lutions and also the embedding free functions 
up to some residual gauge freedom. This matching pro¬ 
cedure simultaneously determines the influence of the 
external gravitational field on the internal structure of 
the neutron star, and also the influence of the neu¬ 
tron star on the external field. The matchings for the 
cases {s,t) = ( 0 , 0 ), ( 0 , 1 ), ( 0 , 2 ), ( 1 , 0 ), ( 1 , 1 ), (U 2 ) have 
been explicitly worked out by Mino et. al. |^, in the 
case where the interior metric describes a black hole. 


In Sec. IV below we will review these matchings and 
also justify some of the assumptions made by Mino et. 
al. concerning the monopole and dipole pieces of the 
fields. We shall need to consider in addition the cases 
(s,f) = (0,3), (0,4), (I, 3), (1,4) and (2,4). 

In order to understanding how the matching works in 
these cases, we need to understand in detail the nature of 
the space of solutions in the internal scheme. We examine 
this issue in the next section. 


III. THE INTERNAL SCHEME : STRUCTURE 
OF THE SOLUTION SPACE 


A. Gauge freedom 

We start by discussing the gauge freedom in the per¬ 
turbation equations. This gauge freedom consists of one 
parameter families of diffeomorphisms 

where V'o is t he id entity map, which act on the 
quantities ( p.2C| ) and (2.21) via the natural pull-back ac¬ 
tion. We can express such a map to O(e^) as 

-Ipe = 22^(1) (£) o 22^(2) (e^) o I?J(3) (e^) o (e^), (3.1) 

where “ for 1 < p < 4 are vector fields and where 
for any vector field t“, I?r(A) : denotes 

the one parameter group of diffeomorphisms generated 
by r“. Since all perturbation quantities vanish at first 
order in e, we can without loss of generality assume that 
^(i)a = 0 . 

Now suppose that S{e) is any one parameter family 
of tensor fields on (we suppress tensor indices), 

which has the expansion 

S{e) = -f -f -f 0{e^). (3.2) 

Then from Eq. & we can calculate the transformation 
properties of the expansion coefficients etc. We 

find 

-f + 0{s^). (3.3) 

where 

§(2) ^ 5'(2) 


and 


gO) ^ 5 ( 3 ) +£^(3)5(0)^ 


5 ( 4 ) = 5 ( 4 ) ^ 

£^(4) ^( 0 ) + ^ 

£ 5 ( 2 ) £^( 2 ) ^( 0 ) 


+£5(2) 5(2). 


(3.4) 

(3.5) 


(3.6) 


Here £ means the Lie derivative. The formulae (3.4) — 
apply when we take S{e) to be any one of the fields 


Pabis), Tab[e), t{e) and t“(e) defined in Sec. HC 1 


We now fix a certain portion of the gauge freedom. We 
can expand the fields t{e) and f“(e) as 

t{e) = -h + 0(£^), (3.7) 


and 


r (e) = t( 0 ) “ + £ 2 ^( 2 ) « + £3^(3) a ^4^(4) a (3 3 ) 

From Eqs. (^.16|) and (|3.4|) - (^) it follows that 


we can 


choose a gauge such that 

i( 2 ) ^ ^(3) ^ ^(4) ^ ^( 2 ) a ^ ^(3) a ^ ^(4) a ^ o_ (3 9 ) 

The residual gauge freedom then consists of the transfor¬ 
mations for which 




(3.10) 


for 2 < p < 4. Thus, in the coordinate system (x\t) of 
Eq. (^), the vector fields are purely spatial, = 0 
and ^^^(x-^,t) = 
































B. The Newtonian case 


To motivate our analyses below, consider first the anal¬ 
ogous Newtonian perturbation theory of a static star . 
The density is given by p = +£‘^p^^+ 

0(£^), the pressure by p = +£^p^^^ +£‘^p^‘^^ + 

0(£^), the velocity by u = -|- -|- 0(e®), and 

the Newtonian potential by $ = -|- 

£4<j)(4) _|_ (9(£5^^ pert urbe d e quatio ns of motion [the 
Newtonian limit of Eqs. ( 2.29|) - ( 2.36| )] are at 0(£^) 


V 2 $( 2 ) _ 47rGp(2) 

Vp(2) + p(2) p(0)vd>(2) 

at 0(£^) 

V2$(3) _ 47rGp(3^ 
Vp(3) + p(3) p(0)vd>(3) 


V • 


pWu(3) 


and at 0(£^) 

V2$(4) _ 47rGp(4) 


V • 


^(o)u(4) 


0 

(3.11) 

0 , 

(3.12) 

0 

(3.13) 

0 

(3.14) 

dt ’ 

(3.15) 

0 

(3.16) 

_p(2)v$(2) 


^ dt 

(3.17) 

9p(3) 

dt ■ 

(3.18) 


The general solution for the perturbation (5$ = $—$(3) 
to the Newtonian potential outside the star is 


= E E 


J—O m— — J ^ 


^ Jm (0 

^ J +1 




Yj^{e,p), (3.19) 


where r, 0, p are the usual spherical polar coordinates. 
We denote the eigenvalue of total angular momentum by 
J here rather than the more conventional I in order to 
facilitate comparison with the relativistic case below. We 
assume that the external gravitational field as parame¬ 
terized by the coefficients bjm varies with £ as 

bjm. = + £^b^Pm + (3.20) 

where each varies over timescales oc l/£, as in the 
relativistic case. Then the coefficients ajm (the star’s 
multipole moments) can be expanded as 

ajm = + £'^af^ -h 0{£^), (3.21) 


where the coefficients are determined from the 
via Eqs. (3.11) - (3.IS). 

There are three well-known properties of the solution 
( |3.19 ) that we will focus on. These properties, suitably 
generalized, carry over to the relativistic analysis and 
will play a crucial role in our matching analysis in Sec. 
IV below. 


• The piece of the solution with total angular mo¬ 
mentum eigenvalue J diverges at large r like r"^. 


• For J > 2, the solutions are completely determined 
by the coefficients at ea ch order in perturba¬ 
tion theory. From Eqs. (3.11) - (3.15) we have 




( 2 ) 

“Sm = kjR 


Jm 

2J+iR3) 


(3.22) 

(3.23) 


where R is the stellar radius and kj is a fixed di¬ 
mensionless constant which is determined by the 
stellar equation of state. At higher orders in £ 
things are a little more complicated. From the form 


of Eqs. (^.16|) - (3.181) it follows that 


consists 


of a piece linear in b^j^ and b^f^, together with a 
term 



J'm' J" m" 


T^Jm 7 ( 2 ) ,( 2 ) 

^ J'm' J" m" ^ J'm' ^ J"m" ’ 


(3.24) 


for some constants Kjj^, j„ . Similarly the solu¬ 
tions inside the star are completely determined by 
the coefficients 6j^. 

• The J = 0,1 pieces of the solution outside the star 
contain physical information only at 0{£'^). New¬ 
tonian conservation of mass forbids solutions with 
0-00 0, and boo is an additive constant to the po¬ 

tential which has no physical consequences. By a 
change of in the origin of coordinates of the from 
r ^ r — ro which makes r = 0 the center of mass 
of the star, one can enforce = 0 for t = 2,3,4. 
This then implies that b[^ = b[^ = 0. However, 
the coefficient 6^^ can be nonzero; this coefficient 
encodes the acceleration of the center of mass of the 
star with respect to inertial reference frames (due to 
couplings of multipole moments of the star to mul¬ 
tipole moments of the external gravitational field 
with J > 2). Similarly, inside the star, J = 0 per¬ 
turbations and J = 1 perturbations to a rotating 
state are driven by interactions with the external 
gravitational field only at order 0(£^). 


C. The relativistic case 


Turn now to the corresponding analysis in the rela¬ 
tivistic case. Let Q denot e t he sp ace of solutions of the 
equations of motion (2.29) - (2.36), consisting of the ten¬ 


sor fields for 2 < t < 4 together with the fluid vari¬ 


ables. The gauge transformations ( ^ ) - (3.6) define an 
equivalence relation on Q. Let S be the set of equivalence 
classes, which is the physical solution space. 

We start by showing that there is a well defined ac¬ 
tion of the rotation group on S. Let i/jr : 
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j^^(NS) isometry of the background spacetime 

corresponding to a rotation R. Under this rotation 


an element 

;( 2 ) 


v'-ah ’ ''‘ab > 

.(3) 


;(4) 


of G gets mapped to 


FromEqs. (|^) - ( |3.q ) 


and using it can be seen that if two el¬ 

ements of G are related by a gauge transformation gener¬ 
ated by the vector fields ^^3) and then their 

images under R are also gauge equivalent, related by the 
gauge transformation generated by 
and Thus the action of the rotation group on 

G extends to an action on S. If Ji is the generator of the 
action of rotations on S, then we can classify elements of 
S in the usual way according to the eigenvalue J(J -I- 1) 
of JiJ’’ and the eigenvalue m oi Jz. 

Next, we define the electric and magnetic pieces of the 
Riemann tensor in the usual way: 


and 


£ab = Racbdf'r 


Bab = ^ecadeR%fi‘'i^ ■ 


(3.25) 


(3.26) 


Here = t“/||t“|| and all quantities are defined with 


respect to the metric ^^nterior) ( 2.1^ ). In our 


chosen gauge (3.9), £ab and Bab will be purely spatial in 
the coordinate system (^^), so we write these tensors as 
£ij and Bij. We will use these variables to discuss the 
large r behavior of the solutions because they are more 
nearly gauge invariant than the metric perturbations. As 
usual we expand these tensors as M 


c _ c(C) 




+ e^£^f+0(6^), (3.27) 


and 


B., = e^B[f + e^Bf + e^B^^ + 0{e^). (3.28) 

These tensors transform via the gauge transformation 
rule (3.4) - ( ^ ) under the spatial gauge transforma¬ 
tions allowed by Eq. ( fl.lO ). The tensors B^j^ and 
are gauge invariant, but the rest are gauge dependent. 

Suppose now that one is given an element of S corre¬ 
sponding to the eigenvalues {J,m). Then it is straight¬ 
forward to show that one can choose a gauge in which 
£ij and Bij are eigenfunctions of J^J* and Jz with the 
corresponding eigenvalues, where Ji is the generator of 
the action of rotations on G- Below, when discussing the 
{J,m) sector of S, we will always assume that such a 
gauge has been chosen. 

To discuss the general form of the solutions it is conve¬ 
nient to use the pure orbital tensor spherical harmonics 
rp\L,jm^Q^ (/?) of Thorne, which are defined in Eqs. (2.40a) 


- (2.40f) of Ref. |41|. Here A, L, J and m are integers 
with A = 0 or 2, 0 < J < oo, — J < m < J, and 


These tensors are complete in the sense that any sym¬ 
metric tensorial function fij{9, ^p) on the unit sphere can 
be expanded as 


fvio,ip)= Y. Y 

A—0,2 J—0m— — J L 


Jm 




XL,Jm 


U 




(3.30) 


As the notation suggests, the tensor T-)) 


XL.Jm f n \ • 

(o', Lp) IS an 

eigenfunction of JiR with eigenvalue J{J + 1) and of Jz 
with eigenvalue m. It is also an eigenfunction of LiU 
with eigenvalue L{L + 1), where Li is the generator of 
the “pure orbital” representation of the rotation group 
on symmetric 2-index tensors |^I|] . Using these tensors 
we can write down the general form for and B^l^ on 
the {J,m) sector of the solution space S: 


U’ = E Uhl 


Jm ij 


(/j) 


A=0,2 L 


and 


E E< 


(i) rpXL^Jm 
Jm ^ ij 




(3.31) 


(3.32) 


A=0.2 L 


Here the coefficients are functions of t and r, where 
(t, r, 0, ip) is the coordinate system ill) , and the sum¬ 
mation over L is over the range ( |3.29 ). Note that the 
eigenvalue L is well defined on G but not on the physical 
solution space S, since the “pure orbital” representation 
of the rotation group does not extend from G to S. 


1. Leading order solutions 


Consider now the general form of the linearized so¬ 
lutions and outside the star. From Eq. 

•'J •'J 


( ^.34 ) these describe a stationary perturbation of the 
Schwarschild spacetime In the flat spacetime limit 
M ^ 0, £\^^ obeys the system of equations 


cC2) _ c(2) * _ Tjic{‘^) _ /-) c(2) _ n 


(3.33) 


Using the methods of Ref. one can show that the 
12 ) 

general solution for is 


7(2) 


= E 


J 

Y - 

J—O m— — J 
J 

Y 


a 


( 2 ) 

Jm rj~t2 J-\-2,Jm 


pJ-\-3 




E 


r(2) J—2 rp2 J—2,Jm 


(3.34) 


J—2m=—J 


J-2<L<J 
1 < L < 3 
L = 2 


for J >2 
for J = 1 
for J = 0 


which is none other than didj 


where is the 
(3.29) O(e^) piece of the Newtonian potential ( ^.19 ). Note that 

the summation in the second term in Eq. (|3.34) starts 
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at J = 2, unlike the Newtonian formula. Also the term 
proportional to is pure gauge. Similarly satisfies 


M2) ^ g(2) * ^ ^ ^(2) = 0, 

[m * y i* j]fc ’ 


(3.35) 


and can be written as 


,( 2 ) 


M2) _ ^Jm rp2J+2,j7n 

2-^ 2-^ rJ+^ b 


J=1 m——J 
OQ J 


E E 

J=2 m— — J 


(2) ^J-2 2^2 J-2,Jm 


Jm 


(3.36) 


In Eq. ( 3.36|) there is no J = 0 ter m in the first sum; such 
a term would satisfy Eqs. ( |3.35 ) but is not present for 
metric perturbations that satisfy the original linearized 
Einstein equation. 

Returning from the flat spacetime limit, the general 

solutions for and in the case M ^ 0 outside the 

( 2 ) ( 2 ) 

star can again be parameterized by coefficients 


.( 2 ) 

-"Jm’ 


and d 


( 2 ) . 
Jm ’ 


Y. Y “S.‘”£r 

J = 0 
J 

(30 J 

*Y Y 


{3.37) 


J—2 m— — J 


and 


^^’ = E E 


^.(2) {\)nJm 
'-'Jm 


J—1 m— — J 
oo J 


+YY 


(3.38} 


J—2 m— — J 


Here and are certain 

fixed linearly independent solutions of the perturbation 
equations in the {J, m) sector whose large r behavior is 
given by @ 


(T) cJm _ 
^ij 


^ rjn2 J-\-2 J JTfl ^ 


M 


jtJ-\-3 ^3 I /pJ-\-A 

T‘f2-2,Jra o{Mr^-^ 




1 




(2 J-\-2 , J TYl 


+ 0 


M 


rpJ-\-3 ^3 I j>J-\-A 

_ j.J-2 ji2 .J-2,Jm 


(3.39) 

(3.40) 

(3.41) 

(3.42) 


These quantities are functions of r, 9 an d (p o nly, a nd ca n 
be expressed as expansions of the form ( 3.31 ) and ( ^.32 ). 
The functions and are gauge independent, 


while a nd are gauge-dependent, their large 

r forms (|3.39 ) and ( R4C|) being achieved only in certain 
gauges In Eq. (|3.37|) we have omitted the J = 1 
term from the first summation since it is pure gauge. 

The above analysis assumes that the (J, m) sec¬ 
tor of the space of stationary linear perturbations off 
Schwarschild, modulo gauge transformations, is of di¬ 
mension 4 for J > 2, of dimension 1 for J = 1, and 
of dimension 1 for J = 0 . Th ese d imensi onali ties are 

strongly suggested by the form ( ^.34 ) and ( 3.36 ) of the 
solutions in the M ^ 0 limit. To rigorously prove these 
dimensionalities in the J > 2 case, one can appeal to the 
Newman-Penrose perturbation formalism. The gauge- 
invariant Newman-Penrose coefficient 'I '4 determines the 
perturbation uniquely up to gauge transformations and 
up to J = 0,1 perturbations In the zero fre¬ 

quency case of interest to us, ^'4 can be written as 

V?) = ^ E -2Y.Jm{0, (fi), (3.43) 

Jm 

where Rjm (r) obeys the zero-frequency Teukolsky equa¬ 
tion 1^. The Teukolsky equation has two complex (or 
four real) linearly independent solutions, since it is a 
second order ordinary differential equation f^ , which 
proves the result. The dimensionalities in the cases 
J = 0,1 can be proved directly; see, e.g.. Ref. f^ . 

The three prope rties of the Newtonian solutions dis¬ 
cussed in Sec. HI B ha ve an alogs in their relativistic coun¬ 
terparts ( 3.37 ) and ( ^3.3^ ). First, in the sector (J,to), 
the Riemann tensor perturbation diverges at large r like 
r'^~^, corresponding to a rate of divergence of the metric 
perturbation (in a suitably chosen gauge) of r'^. Sec¬ 
ond, for J > 2, the solutions are completely determined 
(up to gauge) by the coefficients ^jm which pa¬ 

rameterize the external gravitational perturbations. The 

perturbed fluid equations (|2.29|) and (p.34|) determine the 

^- ' ' - ^f 2 ') ( 2 ) 

response of the star and the coefficients a\' and c\' as 


functions of 6 j 2 j and 






" 4/r2J-l-l j(2) 


-^Jm 

(3.44) 

(3.45) 


This is completely analogous to the Newtonian case 
( |3.22 ) e xcep t that there are now the magnetic type mo- 
ments (3.45|) in addition to the electric type moments 


(|3.44). Here k'j and kj are dimensionless constants de¬ 


termined by the stellar equation of state. Similarly the 
( 2 ) 

fluid variables are determined. 

Finally, consider the status of the J = 0,1 sector s. So - 
lutions of the perturbation equations ( 2.29| ) and (2.34) 
which are purely J = 0 do exist. If we denote the 
background solution with total mass M as 
and Tab[M], then the quantities + e'^SM{t)] 

and Tab[M + e‘^6M{t)] satisfy the perturbation equa¬ 
tions to O(e^) for any choice of 5M{t). However, in 
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this case n o sol uti ons to the high er ord er pe rturbation 
equations ( 2.30| ), ( ^.3l| ) and ( 2.35 ) and ( 2.36 ) exist un¬ 
less 5M{t) = 0. Thus there are no J = 0 perturbations 
at order O(e^). Alternatively, one can argue that J = 0 
perturbations are forbidden at this order in e by conser¬ 
vation of baryon number . 

A similar situation holds in th e J = 1 s ector. The 
J = 1 piece of the solutions (|3.37|) and (|3.38|) is parame- 

( 2 ) (2,') 

terized by the coefficient since the term is pure 
gauge and there are no J = 1 terms in the pieces of the 
solutions that diverge at large r. This coefficient c 


( 2 ) 

Im 


parameterizes a perturbation to a rotating state. The 
second order perturbation equations are satisfied by any 

/n\ 

choice of time dependent coefficient c)^(t), but, as in the 
J = 0 case, the higher order perturbation equations can 
only be satisfied if c^^(t) = 0 |^]. 

To summarize, if we now combine Eqs. (Q, (B , 
( |3.44 ) and (3.45) and revert to using metric perturba¬ 
tions instead of curvature perturbations, we can write 


the solutions for and as 


'^ab — 




J—2 m 


and 


7^(2) _ 

^ab ~ 


oo 

EE[' 

J—2 m 


7 (2) rjiE^J'm 


,(2) 7B,Jm 
^Jm'^ab 


7 ( 2 ) rpB ^JtYI 
^Jm^ab 


(3.46) 


(3.47) 


the O(e^) variables. Hence from Eqs. (3.46) and (3.47) 


we can write 


),(3) _ 0)j(3) 
"'ab ~ "'ab 


rp{^) _ 0rp{3) 

^ab ~ ^ab 


1:1: [i.'f.Uf"'" 


^ Jm'^ab 


J—2 m 
00 




Jm-^ ab 


,(3) f B.Jm 
^Jm"'ab 


7 ( 3 ) rpB^jTfl 

^Jm-^ab 


(3.50) 


(3.51) 


J—2 m 


for some coefficients The quantities and 

here can be regarded as fixed linear functions of 
the time derivatives of hj and or equivalently as 
fixed linear functions of the time derivatives bj^ and 

^Jm- 


Thus, the 0(£^) solutions are completely speci¬ 
fied by giving the coefficients b^jm 

functions of time. 

A similar argument can be us ed at 0{e‘^). Here the 
equations of motion (^.31) and (2.36) for the variables 


^lib contain source terms which are linear in 

time derivatives of the lower order variables, and also 
source terms which are quadratic in the O(e^) variables. 
The argument shows that 


_ of (4) I 

'‘ab ~ '''ab ' 




74 ( 4 ) 07 - 1 ( 4 ) 

^ab ^ab 


J—2 m 
00 


\ ^ \ ^ r^(4) rpE,Jm I ^(4) rpB^Jm 
/ > / > [^Jm-^ab ~^^Jm-^ab 


J—2 m 


(3.52) 
, (3.53) 


Here the electric-type quantities and and 

magnetic-type quantities and are fixed up 

to gauge transformations. 


2. Higher order solutions 

Consider nex t the 0{e^) piece of the solution space S. 
The equations (2.30|) and (2.35) satisfied by the metric 


perturbation and by the 0{e^) fluid variables are 
of the same form as the O(e^) equations, except that 
they contain source terms which are time derivatives of 
the 0{e^) perturbations. Since h'ffj is therefore not a 
vacuum perturb ation of Schwarschild, one cannot use the 
argument of Sec. Ill C 1 to determine the general solution. 
However, we can instead make the following argument. 

r3\ 0 f (3) rp(2) 


Eix a specific 0{e^) solution 
Then, any other solution ) which has 


the same O(e^) part must be of the form 


"■ab — 


0/^(3) 

"ab 


Ah 


(3) 

ab ’ 


_ 07 -’( 3 ) _j_ 

^ab ~ ^ab ^^ab ^ 


(3.48) 

(3.49) 


for some coefficients 5j^, d'j^. 


Jm. Here the quantities 
and depend linearly on b^jl, Jjl, bfl, and 

quadratically on and d^fl ^ as in Eq. ( |3.24 ). 

Note that in Eqs. ( 3.50| ) - ( 3.53 ) we have not specified 




(4) 


yet how the fixed solutions , °Tab , °^ab ^ ^’^d 
are chosen. One has the freedom to add to th ese fun ction s 
any linearized solutions of the form ( |3.46 ) and ( ^.47 ). 
This has the effect of redefining the zero points of the 
variables b^jh, b^jln and It will turn out below 
that for the calculations in this paper we shall not need 
to resolve this ambiguity. 


IV. RESPONSE OF THE STAR TO EXTERNAL 
FIELDS 

In this section we will per form the matching calcula¬ 
tions outlined in Sec. H C 3 above to determine the re¬ 


sponse of the neutron star to the external perturbing 
gravitational fields. 


A. Change in central density 


(3) 

ab 


where from Eqs. ( ^.3C ) and ( ^.35 ) the differences Ah 
and AT^f^ obey the same equations ( 2.29 ) and ( 2.34 ) as 


The change in the structure of the neutron star is de- 

(2) (3) 

scribed by the stress tensor perturbations T^i, and 
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Following Refs. we focus in particular on the 

central density of the star, since changes in the central 
density reflect star-crushing forces that tend to destabi¬ 
lize the star to radial collapse, or “anti-crushing” forces 
that tend to stabilize the star against radial collapse. We 
shall show in this subsection that the change in central 
density depend only on the details of the matchings at 
the orders 0 < t < 2 , and is independent of the details of 
the matchings for t > 3 (see Fig. ^ below). This fact will 
considerably sim plify o ur analysis. 

Following Eq. ( 2.11 ) we expand the density p as 


derivatives of the O(e^) variables and since 

those variables have no J = 0 parts, neither do 
and . Therefore there is no J = 0 part to the O(e^) 
perturbation and hence p^^\t,0) = 0. In addition, the 
quantity Dip^‘^\t,0) appearing in Eq. must vanish, 
as there is no J = 1 piec e to the 0(e^ perturbation. 
Hence, we can rewrite Eq. (4.2) as 


“'( 2 ). 


p = p^^'^ + + O(e^). (4.1) 

We define the central density Pc{t) to be the maximum 
value of the density in the star on a hy persurface of con¬ 
stant t{£). In our chosen gauge (|3.9D , this is the same 
as a hypersurfac e of constant t, where t is the time co¬ 
ordinate in Eq. (2.4). Note that pdt) is not the same as 
p{t,r = 0 ), since the location of the star’s center can be 
changed by the perturbation [|^. Since p^^'^ has a local 
maximum at r = 0 , we find that the change in central 
density is given by 

Spcit) = (t, 0 ) -f s^p^^^ (t, 0 ) 


-he" 


p^^Ht,0)- 


3 D,p^^\t,0)D^p(^\t,0) 
2 D^p(o'){t,0) 


+0{e% 


(4.2) 


Here Di is the spatial derivative associated with the back¬ 
ground metric ■ We will now show that the first two 
terms in the expansion (4.2) are vanishing. 

To prove this result, we modify slightly an argument 
used by Brady and Hughes in a similar context |j^ . The 
change in the stellar structure can be decomposed into 
contributions fro m each of the (J, m) sectors of the per¬ 
turbation [cf. Sec. HIC above]. Therefore we can express 
the density perturbations for 2 < t < 4 as 

p^^\t,r,e,(p) ='^Pjl,it,r)Yjmi0,(p), (4.3) 

Jm 

where (t, r, 9, ip) are Schwarschild-like 

coordinates adapted to the unperturbed, spherical star. 
Since each of the density perturbations p*^*^ is a smooth 
function of position at r = 0 , it follows that p^jl^it, 0 ) = 0 
for all J, m except possibly for the spherically symmetric, 
J = m = 0 sector, as shown by Brady and Hughes [^ . 
Thus, only spherically symmetric sector can change the 
central density. 

Next, from the analysis of Sec. HI C above it follows 
that the O(e^) perturbation contains only J > 2 exci¬ 
tations, and has no J = 0 parts [see Eqs. ( 3.46|) and 


( 3.47 )]. It follows that (t, 0 ) = 0. Similarly, at O(e^), 
the second terms in Eqs. ( t3.50| ) and ( 3.51 ) give a vanish¬ 
ing contribution to p^(t, 0). The terms and 

in Eqs. (3.50|) and (3.51) are linear functions of the time 


Spc{t) = (t, 0) -I- O(e^). 


(4.4) 


Finally, consider the decompositions ( 3.52| ) and ( 3.53| ) 
of the fourth order perturbation variables. The t = A 
matchings can affect only the second terms in these equa¬ 
tions, and those terms will not contribute to the J = 0 
part of the perturbation since they have only J > 2 parts. 
In addition, the only J = 0 piece of the terms and 

are due to the source terms in Eqs. ( 2.33| ) and ( 2.36| ) 
which are quadratic in the O(e^) variables, since the lin¬ 
ear terms will not have any J = 0 part by the argument 
of the last paragraph. Thus, the change in central density 
p^'^^(t, 0 ) is a quadratic function of the O(e^) perturba¬ 
tion variables and is independent of the details of the 
t = 3 and t = 4 matching. 


B. Matching calculations 

Turn now to the matching of the internal and ex¬ 
ternal solutions, which was outlined in Sec. |H C 3 


above. We need to match the metrics ^(^terior) 

order by order in a double expansion m 
the parameters Mjr and r/TZ. The matching scheme is 
illustrated in Fig. 

The matchings for the cases {s,t) = (0,0), (0,1), 
(0, 2), (1, 0), (1, IJ, (1, 2) have been explicitly worked out 
by Mino et. al. p^ ], in the case where the interior metric 
describes a black hole . Their results are also applica¬ 
ble to the neutron star case. They show that vanishes 
identically, and obtain the constraints on the embedding 
free functions that (°’°)(F“(t) describes a geodesic of g^\ 
and that 


(bO)jrO(i) = 0, (bo)jrfe(i) = jfe. 


(4.5) 


Consider now the determination of the metric pertur- 


' (2) L I C2) 

bation . From Sec. pi C| we know that contains 


only J > 2 pieces, and that any (J, m) piece of must 
diverge at large r like r'^. However, from Fig. ^ it can be 

seen that diverges at large r like . Hence, must 

( 2 ,) 

be purely J = 2. There are thus 10 free coefficients 
and 


and ( ^.47[) which determine ■ 
These parameters are determined b y th e (s, t) = (0, 2) 
matching in Fig. which using Eq. (4.5) simply dictates 

f2') 

that the constant asymptotic values of the curvatures fb ’ 
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and associated with agree with the £ij and Btj 

/ D'\ 

of the background metric evaluated on the geodesic. 
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FIG. 1. An illustration of the matching scheme for a single 
body moving in a general vacuum exterior spacetime 5^^^- 
The tensors , h ''^^,... when expanded in power se¬ 

ries in M/r yield the columns in this diagram. Similarly, the 
the tensors g[l\ ^ ^ab "hen expanded in power series in 
r/TZ and when acted on by the pullback map (2.43) yield the 
rows. The matching procedure consists of demanding consis¬ 
tency between the two sets of expansions at each order (s, t) in 
the double expansion, wherein the metric perturbations scale 
as {M/ry {r/TVf. In the external spacetime, the metric per¬ 
turbation is determined by the (s,t) = (1,0) matching, 
and by the (s, t) = (2, 0) and (2,1) matchings. In the 
internal scheme, the metric perturbation vanishes identi¬ 
cally, is determined by the (s, t) = (0, 2) matching, 
is determined by the (0,3) and (1,3) matchings, and by 
the (0,4), (1,4) and (2,4) matchings. 


Turn next to the perturbations and . The 
free parameters in these perturbations are the quantities 
and in Eqs. (|^ - (|T^) , and are 
deter mined by the matching scheme. As explained in Sec. 
iV A above, the values of these parameters will not affect 
the change in central density of the neutron star and 
thus are unimportant for our purposes. For completeness 
we briefly mention how these parameters are determined. 


First, since iSp diverges like at large r, it contains only 
J = 3 and J = 2 pieces. The J = 3 piece is determined 
by the (s,t) = (0,3) matching, and depends linearly on 
the spatial derivatives of the electric and magnetic parts 
of the Weyl tensor of the background spacetime gp\ 
evaluated on the geodesic. The J = 2 piece is determined 
by the (1,3) matching, and depends on the first order 
perturbation fipj in the external spacetime as well as 


on the background metric ■ Specifically, hP is first 
determined by the (1,0) matching |^, and from this 
one can calculate the (1, 3) element of the matrix in Fig. 
1^, and hence infer the J = 2 piece of hP. Note that 
hP thus depends on the geometry of gp^ not just in a 
neighborhood of the geodesic, but also non-locally Q. 
In a similar way hP contains pieces with 0 < J < 4 in 


the term °hP in Fq. (|3.52|) which are independent of the 
t = 4 matchin gs, a nd pieces with 2 < J < 4 in the second 
term in Fq. ( p.52| ) which are determined by the t = 4 
matchings. These additional 2 < J < 4 pieces depend on 
hP and hP in addition to gp^ and are determined by 
the (0,4), (1,4) and (2,4) matchings. 

Returning to the 0{e^) perturbation, it follows from 
the arguments of Sec. IV A that the leading order change 
(|4.4D in central density depends quadratically on , 
and hence quadratically on £ij and Bij , the curvatures of 


/ D\ 

the external background spacetime g\^^ evaluated on the 
worldline. Furthermore, since all the relevant equations 
are elliptic, the dependence of 6pc on and Bij is local 
in time. Hence, invariance under rotations yields 


= ci£„(t)£:*^(t) + C2B,yt)B^^{t) + 

Pc 

(4.6) 


where ci and C 2 are constants of dimension (length)”^ 
that depend on M, R and on the equation of state. (A 
cross term between the £ij and Bij fields is forbidden by 
parity argument s.) 

Thus, Fqs. ( ^.29| ) - ( 2.31 ) have the same structure 
as their N ewton ia n an d post-Newtonian count er parts . 
Equations ( 2.29| ), ( ^.3C| ) together with Eqs. ( 2.34 ), ( 2.35| ) 
describe distortions of the star induced by the external 
tidal fields, where there is no change in the star’s cen¬ 
tral density and n o ex citation of the star’s J = m = 0 
modes. Equation (2.31) describes the second-order effect 
of the leading-order tidal field on the star’s struc¬ 
ture. Exactly as in Newtonian gravity, it is this second- 
order effect of the leading-order tidal field that excites 
the spherically symmetric, radial modes of the star and 
changes the star’s central density and angle-averaged ra¬ 
dius. 

It is possible to understand this result in a fairly sim¬ 
ple, intuitive way. Consider first a weakly self-gravitating 
body in an external gravitational field in the tidal limit 
R<^TZ. Simply analyzing the d ynam ics of the test body 
using Fermi-normal coordinates {2^) allows one to imme¬ 
diately conclude that the effect of the external field on the 
body’s internal dynamics must scale as 7?.“^; this is just 
the equivalence principle. The fact that there is no spher¬ 
ically symmetric interaction at this order (that is, that 
the interaction can be called a “tidal” interaction) follows 
from algebraic properties of general relativity - the rela¬ 
tivistic generalization of the familiar fact that the trace 
of the Newtonian tidal force tensor d^^jdx'^dx^ vanishes. 
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As a consequence of this vanishing of the spherically sym¬ 
metric interaction, all radial crushing or anti-crushing 
forces must scale as At first sight the above argu¬ 

ment does not apply to a strongly self-gravitating body. 
However, the essence of the argument can be carried 
through. What is relevant for determining the interac¬ 
tion between the body and the external field are Ein¬ 
stein’s equations in the matching region i? ^ r <C 7?. (the 
body’s “local asymptotic rest frame” ||^). The asymp¬ 
totic value of the external spacetime’s curvature tensors 
in this region (which is the region r ^ 0 as seen in the 
external spacetime) act as a source for the interaction, 
and their scaling (oc TZ~'^) and algebraic properties de¬ 
termine the nature of the interaction in the same way as 
for a weakly self-gravitating body. 

To conclude, we have shown that, just as in Newtonian 
gravity, the leading order change in the central density 
of a fully relativistic spherical star freely falling in an 
external vacuum gravitational field is given by the star’s 
second-order response to the leading order, external tidal 
field. 


V. MODIFIED MATCHED ASYMPTOTIC 
EXPANSION METHOD APPLICABLE TO 
NEUTRON STARS IN A BINARY 


The analysis of Secs. i-0 assumes that the space- 
time outside the body is vacuum. Thi s assumption en¬ 
tered in the equation of motion ( ^.38 ) satisfied by the 
metric perturbation h^^fj , which is used to derive the fact 
that the star travels along a geodesic of the background 
metric g^^^ to leading order [cf. Eq. (^A^ above]. It is 


possible to modify the analysis of Secs, y - to ac¬ 
commodate two freely falling bodies, for examples two 
neutron stars moving in the vicinity of a supermassive 
black hole. In this case the metric perturbation h^^^j de¬ 
scribes the linearized gravitational interactions of the two 
neutron stars, and one must solve simultaneously for the 
external metric perturbations, for two sets of internal 
metric perturbations, one for each star, and likewise for 
two sets of embedding functions. The scheme allows one 
to derive, for example, the equations of motion of two 
“point particles” interacting via their linearized gravita¬ 
tional fields. 

Such a calculational scheme is applicable in principle 
to an isolated neutron star binary, but is poorly adapted 

/ D\ 

to that situation. Since the background metric ' is 
a flat Minkowski metric in this context, the calculations 
of Secs. Q - [V of the leading order change in central 
density are not applicable. Moreover, to achieve our goal 
of determining the scaling of the change in central den¬ 
sity with the parameters e and a discussed in Sec. [A, 
one should describe the gravitational interactions of the 
two neutron stars not by metric perturbations h^^jj and 
^ab ’ rather in terms of a post-Newtonian expansion. 


Eor these reasons, in this section we outline a modified, 
matched asymptotic expansion calculational method in 
which metric perturbations in an internal scheme are 
matched onto post-Newtonian quantities in an external 
scheme. The modified method will allow us to calculate 
the change in central density for neutron stars in a binary. 

As in Sec. |i| above, the method consists of an inter¬ 
nal scheme, an external scheme, and a matching scheme. 
A key difference is that there are two sets of internal 
schemes and two sets of matchings, one for each star, all 
of which must be solved self-consistently. 


A. The external scheme 


Let be the external manifold in which the neu¬ 

tron stars move. The gravitational field is described in 
the external scheme by the standard post-Newtonian ex¬ 
pansion in vacuum. Thus the zeroth order, background 
solution is just a Newtonian spacetime, instead of the 

/ D\ 

vacuum Lorentzian metric ' we had previously. The 


metric perturbations h^fl are replaced by post-Newtonian 
fields, post-post-Newtonian fields, etc Q. 

Now, any Newtonian spacetime can be characterized 
by a lengthscale C and a massscale M, such that the 
typical value of the Newtonian potential is ^ M/C and 
such that the local radius of curvature TZc is given by 
TZ~^ ^ M/C^. In our example of a neutron star binary, 
C will be just the orbital separation L. The post-s/2- 
Newtonian fields, for s = 0,1,2..., scale as but 

have no definite scaling with respect to C. T his is anal¬ 
ogous to the behavior of the expansion ( 2.37 ). As is well 
known, the post-Newtonian fields corresponding to odd 
values of s vanish identically for s = 1, 3 and start at 
s = 5. 


In a suitable coordinate system a;“ = (t, x*), the met¬ 
ric up to post-l-Newtonian order can be written in the 
standard form 


ds^ = ~ [l + 2e^<i)(a;*, et) et)'^ 2e^d'(a;*, et) 

-HO(e®)] dP + 2dx^dt [e^Wi{x\et) + 0(£®)] 
-\-dx’'dx^ [Sij — 2e^^{x’',et)Sij 0{£‘^)\ , (5.1) 


where £ oc ^/M is a formal expansion parameter. Equiv¬ 
alently, by making a gauge change t —> t/e, the metric 
can be written as 


ds = -j- 

£^ L 


l + 2e'^^{x\t) + 2e'^^{x\t)'^+ 2 £^-^{x\T) 

-HO(£®)j dP + 2dx^di [e^Wi{x\t) -t 0(£^)] 
+dx^dx^ [Sij - 2£‘^<^{x\t)S^J + 0 (£^)] . ( 5 . 2 ) 


In vacuum one can pick a gauge in which the potentials 
$, and Wi obey the equations 

-h $ = 0, (5.3) 
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where is the Laplacian of flat space. 

The post-Newtonian expansion can also be described 
in a coordinate-free way |^|. Let 5 “*'(e) be a one pa¬ 


rameter family of vacuum metrics which are C'^ in 
in a neighborhood of £ = 0, such that the limit £ —> 0 
of fl°^(£) exists and is of signature (0, -|-, -I-, -|-). [The for¬ 
mula (pj) gives an approximate version of such a family]. 
Then there exist tensor fields h°'^, ta and a connection Da 
such that 


( 


g-\e) = -F 0(/) 
s'^gabis) = -tah + 0(£^) 
Vaie)=Da + Oie^). 


(5.4) 

(5.5) 

(5.6) 


The quantities Da, h°‘^ and ta comprise a Newtonian 
spacetime. In Newtonian or “inertial” coordinate sys¬ 
tems (a:*,t), these quantities are given by ta = {dt)a, 
= 5ij{d/dx'‘Y{d/dx^)^, and Da is given in terms 
of the Newtonian potential $ by the only non-vanishing 
connection coefficient being The Newtonian 

fields satisfy the additional relations 


Dah = Dah^^ = = 0 , 


(5.7) 


where R\ad i® curvature of the connection D 


The 

last of the relations (5.7) is just the the limit £ ^ 0 of the 
identity ~ higher order correc¬ 


tion terms in Eqs. (^) - (^) of order 0(£®) collectively 
describe the post-s/2-Newtonian fields; we shall not need 
the precise forms of these fields here gj] . 


B. The internal scheme 


The intern al scheme is very similar to that described 
in Sec. IIC 1 above. The interior metric still is given by 
Eq. (^.^), however now each term £‘h^2 expan¬ 
sion scales as rather than 77“*, and has no definite 

scaling with respect to M. The reason that we need to in¬ 
clude half-integral powers of l/£ is that time derivatives 
will scale like 

The equations of motion ( 2.291 ) - ( 2.36| ) are modified by 
the following two considerations. First, time-derivatives 
now scale as £^ rather than e, which modifies the form 
of the perturbation analysis. Hence, a time derivative of 
the held will enter as a source term in the equation 


^ab 

for the field h 


(i+3) 

ab 


The perturbations with t odd are 


due entirely to such time-derivative source terms. Since it 
turns out that the first non-vanishing perturbation occurs 
at 0(£®), the first non-vanishing with t odd occurs 

for t = 9. The second consideration is that we need now 
to consider the perturbations up to order t = 12 instead 
of up to order t = 4. This is not as complex as it seems 
since it suffices to consider even values of t until t = 9. 
In addition, as already mentioned the first non-vanishing 


perturbation is 


( 6 ) 

ab 


The structure of the resulting equa tio ns of motion is 
closely analogous to that of Eqs. ( ^.29|) - ( |2.3(j ), and the 
gener al solut ions fo llow the same pattern as the solutions 
(p.46|) , ( tl.47 ) and ( 3.5C ) - ( |3.53|) . As before, the crucial 
aspect of the solutions that we shall use is that the first 
non-vanishing perturbation in the J = 0 sector cannot 
arise directly from matching to the external scheme, but 
must arise from a source term that is quadratic in a lower 
order held. 

In what follows, we neglect entirely odd values of t. 
This is justified since, just as in Sec. IV B above, all time- 


derivative-generated terms have no J = 0 parts and will 
not contribute to the leading order change in central den¬ 
sity. 


C. The matching scheme 

The constr uction of the matching scheme parallels that 
given in Sec. IIC 3| above. One needs to specify an em¬ 
bedding 






(5.8) 


of the neutron star spacetime into the external spacetime, 
i.e., a mapping between the asymptotically Lorentzian 
coordinates a;“ of and the coordinates x°‘ of 


We can write this mapping as [compare Eq. (2.43) above] 


x°‘ = x°‘{t,xY 

p=0 ^ r—0 Ir 


(5.9) 


where A denotes the mu lti-in dex (ii,..., v) and x^^ = 
x^^x’'^ ... x^''. Equation ( ^.43| ) is a Taylor expansion of 
Lp in terms both of the spatial coordinates x® at each 
fixed t, and also in terms of the parameter y^M/C. The 
terms with p odd are vanishing for p = 1,3; the first non¬ 
vanishing terms with p odd start at p = 5. Thus, the 
worldline of the center of the neutron star gets mapped 
onto the worldline 


x“(i) = (0'0)^“(t). 


M 

T 


+ 0[{M/CY]. (5.10) 


The matching procedure described below can be used 
to show that the first term in Eq. ( 5.10| ) represents a 
worldline satisfying Newtonian equations of motion, that 
the second term is the first post-Newtonian, point mass 
correction, etc. 

Next, from the interior metric qj^ 

one construct the metric 


/ -i\ (J 

gab=[p )^x*9a. 


(INTERIOR) 


(5.11) 


on Here v: ^ is the mapping de¬ 

fined after Eq. ( 2.17 ) above, with the parameter e chosen 
to have the value V^M. The metric will 
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have a dependence on '/M like that of the metric ( |5.2| ) 
has on e. Now perform a double power series expansion 
of the metric gab, its inverse 5 “^ and its associated con¬ 
nection Vq in t erms of the parameters M/r and r/£, 
as in Eq. ( 2.40| ) above. These double power series ex¬ 
pansion must be consistent, order by order, with expan¬ 
sions in powers of r/C of the post-s/2-Newtonian fields 
for s = 0, 2,4, 5,... of the external scheme. As before, 
demanding such consistency determines both the embed¬ 
ding free functions and the appropriate solutions in the 
internal and external schemes, up to some gauge freedom. 
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FIG. 2. An illustration of the matching scheme for a body 
moving in a spacetime described by a post-Newtonian approx¬ 
imation. In the external spacetime, the post-s/2-Newtonian 
fields for s = 0,2,4, when expanded in power series in r/C, 
yield the rows in this diagram. The sum of 

• when acted on by the pullback map ( |5.1l| ) yields a 
tensor field gab{M) on the external spacetime. When one con¬ 
siders the triple (Mga6(M),p“*’(M), Va(M)), and expands 
this triple order by order in M, one obtains at the 

post-s/2-Newtonian fields, for s = 0,2,4, 5.... The expan¬ 
sions in powers of r/C of these fields must agree with those of 
the external spacetime post-Newtonian fields. In the internal 
scheme, the metric perturbations and h!'^^ vanish identi¬ 
cally, and is determined by the (s,t) = (0,6) matching. 
The first change in central density occurs in . 


Consider for example the case s = 0, which is just the 
matching onto a Newtonian spacetime. The limit M —> 
0 of the fields (jab, 5 “^ and Va yields yields Newtonian 
fields ia and Da, cf. Eqs. (^) - (^) above. The 
expansion in powers of r/C of these fields must agree, 
order by order, with the expansions in powers of r/C 
of the Newtonian fields ta, and Da of t he e xternal 
scheme. If we choose the coordinate system (2.4) in the 
internal scheme and a Newtonian or inertial coordinate 
system (x*,t) in the internal scheme, then the fields h°'^, 
and ta, ia will coincide to all orders in r/C if we 


choose the embedding ( |5.8| ) to be of the form i{x^ ,t) = t, 
x*(x£t) = X® -I- /*(t). Then one just has to match the 
Newtonian potentials order by order in r/C. 

This matching to zeroth order in r/C [i.e., the (s,t) = 
(0, 0) matching] yields the usual Newtonian point-mass 
equations of motion. Hence, when one solves for the em¬ 
bedding functions of both stars and for the interior and 
exterior gravitational fields up to order {s,t) = (0, 0), the 
result is two static spherical neutron stars moving along 
their Newtonian orbits. This then serves as the starting 
point for calculating higher order perturbations. 

This matching procedure is illustrated schematically in 
Fig. 1^. Our arguments below will be independent of the 
details of the post-Newtonian and post-post-Newtonian 
matchings, so we do not need to describe those here. 


D. Change in central density 


The leading order change in central density Spc of a 
neutron star in a binary can now be dedu ced by argu¬ 
ments analogous to those given in Sec. IV B above. First, 


in the internal scheme w e nee d consider only even values 
of t, as explained in Sec. V_B. Next, consider the pertur¬ 


bations and h'j/'f/. If these perturbations were to be 
non-zero (or not pure gauge), then they would need to di¬ 
verge at large r like for some J > 2, from the analysis 
in Sec. HI. But from Fig. ||, ^ const and oc r 

at large r. Hence the first non-vanishing perturbation is 
h^ab ^ which is purely J = 2. Note that this perturbation 
is determined entirely by the exterior Newtonian fields 
from the {s,t) = (0,6) matching. 

Since the first non-vanishing perturbation is the 


;(4) 


first change in central density must scale as (/i^^ 
the same arguments as before, so we must have 

Spc 1 
-cx —^ 

Pc 


(6)'i2 


by 


(5.12) 


at large C. The constant of proportionality in Eq. ( 5.12| ) 
must have dimensions of (length)®, and the only relevant 
dimensionful parameters are the mass M of the star and 
its radius R. We can therefore write 


^-^ = F[M/R](^ =F(e)a® + 0(a"), (5.13) 


where F is a dimensionless function of a dimensionless 
variable, and we have used the dimensionless variables 
a = R/L and e = M/R discussed in the Introduction. 
The function F will depend on the equation of state. 


VI. DISCUSSION AND CONCLUSIONS 

As discussed in the Introduction, the result (HI) is 
in disagreement with Refs. ||^,|| in the regime a <C 1 
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where there should be agreement, since the scaling found 
in the numerical simulations is 5pc/Pc oc as a ^ 0. 
This disagreement in scaling at large orbital separations 
is strong evidence that the star crushing effect seen in 
Refs. [110 is not physical. 

In addition, at the location of the instability seen by 
Wilson et. al. , the dimensionless parameter a® is of orde r 
10“^, so one might expect the perturbative result ( 5.13| ) 
to be a good approximation. For F (e) of order unity, the 
fractional change in central density ( 5.13 ) at this location 
is ^ 1000 times smaller than that seen in the numerical 
simulations 

The analysis of this paper does not determine the sign 
of the function F{e). However, we can make the following 
argument. Let us expand this function as a power series, 
which yields an expression of the form 


he 


= Foa^ [1 + Fie + + ...] 


( 6 . 1 ) 


for some constant coefficients Fb, Fi, F 2 etc. Now, it is 
known from Lai’s Newtonian analysis that the coefficient 
Fo is negative |@,^. Therefore, unless the dimensionless 
coefficients Fi, F 2 etc are negative and large (~ —7 for 
Fi), which seems unlikely, the change in central density 
due to the tidal interaction will be negative. Thus, the 
stars angle-averaged radius will increase and the star will 
be more stable and not less stable to radial collapse. 

Note that our analysis assumes that the neutron stars 
are not spinning in their local rest frames, although the 
dragging of inertial frames means that they will spin 
slightly with respect to distant stars; this effect is in- 
corp orat ed in the choice of Fermi-normal coordinates in 
Eq. dJ). 


Note also that our analysis assumes that the external 
tidal fields are slowly varying compared to the dynam¬ 
ical timescales of the neutron star, so that the stellar 
modes adiabatically follow their driving forces. This fol¬ 
lows from the assumption a <C 1 underlying the pertur¬ 
bation expansion. Clearly this adiabatic approximation 
breaks down at a ~ 1. However, one can ask how well we 
might expect the adiabatic approximation to be working 
at the location of the instability. For the neutron star 
/-modes and p-modes, the idealization is a good approx¬ 
imation: From the point of view of the companion star, 
the timescale over which the external tidal fields is vary¬ 
ing is T ~ hL^/M. This timescale is long compared to 
the internal dynamical time of the neutron star (the char¬ 


acteristic timescale of the /-modes) ~ i/IFJm. Hence, 
these modes will equilibrate rapidly in response to the 
tidal perturbations and the approximation is fairly good. 
The p-modes of the neutron star on the other hand have 
frequencies <C and are resonantly excited [ po|J^ 

in the Newtonian approximation near the end of the in¬ 
spiral; for these modes the approximation is inappropri¬ 
ate. This is a limitation of the applicability of our anal¬ 
ysis. However, it seems unlikely that resonant p-mode 
excitations could be responsible for the star crushing ef¬ 
fect seen in the Wilson-Mathews-Marronetti simulations, 


since the predicted amplitudes of excitation in the New¬ 
tonian approximation are fairly small . 

Finally, consider the implications of our result for the 
gravitational wave signature of the inspiral. Previous 
analyses [ pfljT^ jTlJl^ have shown that to Newtonian or¬ 
der, the effect of the finite size of the neutron stars on 
the gravitational waveform is small, and in fact is negli¬ 
gible for the purposes of signal detection. For the initial 
LIGO interferometers, ~ 95% of the signal-to-noise ra¬ 
tio will have been accumulated before the gravitational 
wave frequency reaches 400 Hz [see, e.g., Eq. (2.23) of 
Ref. [^]. The total accumulated phase error in the sig¬ 
nal predicted by Newtonian tidal interactions by 400 Hz 
is <C 1; see Eq. (7.5) of Ref. |Q| and Ref. [|^. We have 
shown using a fully relativistic treatment of tidal inter¬ 
actions that, as one would expect, the Newtonian predic¬ 
tions for tidal interactions are valid except for correction 
terms of order M/R and M/L, that is, fractional correc¬ 
tions < 1. Therefore, the effects of tidal interactions and 
of the finite size of the neutron stars will be unimportant 
for signal detection. 

As this paper was being completed, we learned of a 
similar analysis by Thorne [Q. Thorne shows that the 
change in central density is oc a®, but also that it is 
always negative, stabilizing the neutron star. 
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TABLE I. A summary of the orders at which various physical effects occur. The parameter a is the tidal expansion parameter 
(stellar radius)/(orbital separation), and e is the strength of internal gravity in each star, e =(stellar mass)/(stellar radius). A 
perturbation to the orbital motion is said to scale as if the ratio of the perturbing force to the Newtonian orbital force 

is of this order. Similarly, the distortion of each star is said to be of order if the dimensionless measure of distortion 

{Ixx — Iyv)/Ixx is of this order. Note that to post-l-Newtonian order, all the physical effects scale the same way with a as 
their Newtonian counterparts. The pattern for the stellar modes is repeated at all higher orders in e. The entries in the orbital 
motion columns can be deduced from the entries in the stellar modes columns; if 5£ is the energy transferred between the stellar 
modes and the orbital motion, then the fractional change in the orbital motion is ~ LSS/M'^, while the fractional distortion of 
the star is ~ R5£/M^, one power of a smaller. 
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